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q . Abstract. 

We provide geometric quantization of a completely integrable Hamiltonian system in the action- 
angle variables around an invariant torus with respect to the angle polarization. The carrier 
^ | space of this quantization is the pre-Hilbert space of smooth complex functions on the torus. A 
fvq ' Hamiltonian of a completely integrable system in this carrier space has a countable spectrum. If 
it is degenerate, its eigenvalues are countably degenerate. We study nonadiabatic perturbations 
of this Hamiltonian by a term depending on classical time-dependent parameters. It is originated 
by a connection on the parameter space, and is linear in the temporal derivatives of parameters. 
One can choose it commuting with a degenerate Hamiltonian of a completely integrable system. 
Then the corresponding evolution operator acts in the eigenspaces of this Hamiltonian, and is an 
Q ' operator of parallel displacement along a curve in a parameter space. 
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1 Introduction 

By the well-known KAM theorem [|T], [OJ], there are perturbations of a Hamiltonian of a 



completely integrable Hamiltonian system around an invariant torus such that trajectories 
q-i of a perturbed Hamiltonian system do not leave a compact neighbourhood of this torus. 

Here, we show that one can obtain the similar effect by means of a suitable control operator, 
though trajectories of the control system need not live on tori. The perturbed system that 
we construct depends on parameters. A Hamiltonian of such a system is a sum of a dynamic 
Hamiltonian of the initial completely integrable system and a perturbation term which is 
linear in the temporal derivatives of parameters. The key point is that integration of this 
term over time through a parameter function depends only on a trajectory of this function in 
a parameter space, but not on its parameterization by time. The corresponding evolution 
operator can be treated as parallel displacement with respect to some connection on a 
parameter space, and plays the role of a classical geometric Berry factor. One can use it 
as a control operator in a perturbed completely integrable Hamiltonian system. The goal 
is quantization of such a system. 

Note that, at present, geometric factor phenomena in quantum systems with classical 
parameters attract special attention in connection with holonomic quantum computation 
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T7|, E3J. This approach to quantum computing is based on the generalization of Berry's 



adiabatic phase to the non-Abelian case corresponding to adiabatically driving an n-fold 
degenerate eigenstate of a Hamiltonian over the parameter manifold ||21|| . Information is 
encoded in this degenerate state, while the parameter manifold plays the role of a control 
parameter space. 

There are different approaches to quantization of autonomous completely integrable 
Hamiltonian systems || [J. An advantage of geometric quantization is that it remains 
equivalent under symplectic isomorphisms. One has studied geometric quantization of 
an autonomous completely integrable Hamiltonian system around an invariant torus with 
respect to polarization generated by Hamiltonian vector fields of first integrals of this system 



14] , [15i . The problem is that the associated quantum algebra contains functions which are 



not globally defined, and that elements of its carrier space fail to be smooth sections of the 
quantum bundle. We choose a different polarization. 

Let (Z, Q) be a 2m- dimensional symplectic manifold which admits m smooth functions 
Fk, k = 1, . . . , m, which are pairwise in involution and whose differentials dFk are linearly 
independent almost everywhere. This is called a completely integrable Hamiltonian system. 
Let M be its (connected) compact manifold, i.e., all Ff, are constant on M. By the well- 
known theorem [|I], [LT| , if dFk 7^ on M, there exists a small neighbourhood of M which is 
isomorphic to the symplectic annulus 

W = V x T m , (1) 

where V C M m is a nonempty (open, contractible) domain and T m is an m-dimensional 
torus. The W is equipped with the action-angle coordinates (Ik, 4> k mod 2ir). With respect 
to these coordinates, the symplectic form on W reads 

n = dI k A d(f) k , (2) 

and all Fk are functions of action coordinates (Ik) only. A Hamiltonian on W is an arbitrary 
analytic function TC of action coordinates Ik- The corresponding Hamilton equation reads 

i k = 0, <p k = d k n. (3) 

In order to geometrically quantize the symplectic manifold (W, Q), we choose the angle 
polarization spanned by the almost-Hamiltonian vector fields d k of angle variables |J. The 
associated quantum algebra A consists of functions which are affine in action variables I k . 
We obtain the continuum set of its nonequivalent representations by first order differential 
operators in the separable pre-Hilbert space C°°(T m ) of smooth complex functions on the 
torus T m . This set is indexed by homomorphisms of the de Rham cohomology group 
H l (T m ) of T m to the circle group U(l), i.e., by collections [A&], k = l,...,m, of real 
numbers Xk G [0, 1). In particular, the action operators read 

h = -idk - A fc . (4) 
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By virtue of the multidimensional Fourier theorem M, the orthonormal basis for C°°(T m ) 
consists of functions 

VV) = exp[i(n r r )], (n r ) = (n x , ...,n m )e Z m . (5) 
These are eigenvectors 

of the action operators (0J) and, consequently, of a Hamiltonian 

?#(n r ) = W(%-)VV) = - A i)^K)- (6) 

In particular, if a Hamiltonian of a completely integrable system is independent on some 
action variables, its eigenstates are countably degenerate. In the framework of holonomic 
quantum computation, such a degenerate state can be utilized for encoding information. 
The goal is to build a holonomic control operator acting in this state. We construct it 
depending on classical time-dependent parameters. 



2 Quantization of a completely integrable Hamiltonian system 



We follow the standard geometric quantization procedure |2], |20], p2| . Since the symplectic 
form Q (|2]) is exact, the prequantum bundle is a trivial complex line bundle C — ► W. Let 
its trivialization 

C = W x C (7) 

hold fixed. Any other trivialization leads to equivalent quantization of W. Given the 
associated bundle coordinates (Ik,4> k ,c), c G C, on C (0), one can treat its sections as 
smooth complex functions on W. 

The Konstant-Souriau prequantization formula associates to each smooth real function 
/ G C°°(W) on W the first order differential operator 

f = -iV* t +f (8) 

on sections of C. Here 

$ f = d k fd k -d k fd k 

is the Hamiltonian vector field of /, and V is the covariant differential with respect to a 
[/(l)-principal connection AonC whose curvature form obeys the prequantization condition 
R = iQ. Such a connection reads 

A = A + icl k d(f) k <g> d c , (9) 
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where A is a flat [/(l)-principal connection on C. The classes of gauge conjugated flat 
principal connections on C are indexed by the set M m /Z m of homomorphisms of the de 
Rham cohomology group 



H\W) = H\T m ) =R m 
of the annulus W (0]) to Z7 (1) 0. Let us choose the representatives 

A [X k ] = dl k <g> d k + d(j) k ® (d fe + zA fc c<9 c ), A fc e [0, 1), 
of these classes. Then the connection A (0) up to gauge conjugation reads 

A[A fc ] = dl k ®d k + d<\> k ® (d k + i(I k + X k )cd c ). (10) 



For the sake of simplicity, let Afc in the expression (|10| ) be arbitrary real numbers, but we 
will bear in mind that connections A[X k ] and A[\' k ] with A^ — A' fe G Z are gauge conjugated. 
Given a connection (|10|), the prequantization operators (H) read 



f = -i# f + (f-(I k + \ k )&f). (11) 

Since the divergence of any Hamiltonian vector field with respect to canonical coordinates 
vanishes, the prequantization operators / also keep their form ( |TTD on sections of the quan- 
tum bundle C(g)X>i/2, where Vi/ 2 — > W is a metalinear bundle, whose sections are half-forms 
on 

Let us choose the manifested angle polarization. It is the vertical tangent bundle Vn of 
the fiber bundle 

7T : V x T m -> T m 

spanned by the vectors 9 fc . The corresponding quantum algebra .A C C°°(W) consists of 
affine functions 

f = a k {4r)h + b{4T) (12) 



of action coordinates I k . The carrier space S of its representation ( |TT| ) consists of sections 
p of the quantum bundle C ® T>i/ 2 ^ W oi compact support which obey the condition 
V,?p = for any Hamiltonian vector field d subordinate to the polarization Vn. This 
condition reads 

d k fd k p = 0, VfeC°°{W). 

It follows that elements of £ are independent of action variables and, consequently, fail to 
be of compact support, unless p = 0, i.e., E reduces to p = 0. 
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Therefore, let us modify the standard quantization procedure as follows [0. Given an 
imbedding 

i T : T m -> V x T m , 

let Ct = i^C be the pull-back of the prequantum bundle C (|7|) onto the torus T m . It is a 
trivial complex line bundle Ct = T m x C whose sections are smooth complex functions on 
T m . Let 

A[X k ] = i T A[\ k ] = d<f) k <g> (d fe + + A fc )c<9 c ) 

be the pull-back of the connection A[X k ] ( |I~0D onto Or, and let V denote the corresponding 
covariant differential. Let T>t be a metalinear bundle of complex half-forms on the torus 
T m . It admits the canonical lift of any vector field r on T m , and the corresponding Lie 
derivative of its sections reads 

L r = r k d k + X -d k r k . 
Let us consider the tensor product 

Y = C T ® V T -f T m . (13) 
Since the Hamiltonian vector fields 

# f = a k d k -(I r d k a r + d k b)d k 

of functions / (fOP are projectable onto T m , one can associate to each f <E A the first order 
differential operator 

/ = (-iV^ + /) ® Id + Id ® L 7r ^ / = -m fe a fc - l -d k a k - a k \ k + 6 (14) 

on sections of Y . A direct computation shows that the operators (|14[) obey the Dirac 
condition 

[/,/'] = -^{/7/ / }. 

Sections s of K — >■ T m constitute a pre-Hilbert space St with respect to the nondegenerate 
Hermitian form 

(s\s')= (J^j J ss\ s,s' eS T . 

Then / are Hermitian operators in They provide the desired geometric quantization 
of a completely integrable Hamiltonian system on the annulus W (|l|). 



5 



Of course, this quantization depends on the choice of a connection A[Xk] flTU| ) and a 
metalinear bundle T*t-, which need not be trivial. 

If T>t is trivial, sections of the quantum bundle Y — > T m (|13|) obey the transformation 
rule 

+ 2vr) = s{<p k ) 

for all indices k. They are naturally complex smooth functions on T m , i.e., the carrier space 
St coincides with the pre-Hilbert space C°°(T m ) of smooth complex functions on T m . Its 
orthonormal basis consists of the functions ^(rv) ©■ The ac ti° n operators Ik ( |l4l) in this 
space take the form (f|). Other elements of the algebra ^4 are decomposed into the pull-back 
functions 7r*ip(n r ) on W which act on C°°(T m ) by multiplications 

^Vv^K) = ^(nrjV'K) = Vv+O- ( 15 ) 



If XV is a nontrivial metalinear bundle, sections of the quantum bundle Y — > T m ( [13D 
obey the transformation rule 

p T (^ + 2tt) = -pt(^) (16) 

for some indices j. In this case, the orthonormal basis for the pre-Hilbert space St can be 
represented by double-valued complex functions 

VW.) = exp[iOi0 J + (rij + (17) 
on T m . They are eigenvectors 

of the operators while the pull-back functions ir*ip(n r ) act on the basis (0) by the 

above law ([15|) . It follows that the representation of the quantum algebra A determined 
by the connection A[\k] ( |i0|) in the space of sections ([T6| ) of a nontrivial quantum bundle 
K ( |13|) is equivalent to its representation determined by the connection A[Xi, Xj — |] in the 
space C°°(T m ) of smooth complex functions on T m . Therefore, one can restrict the study 
of representations of the quantum algebra A to its representations in C°°(T m ) defined by 
different connections (|10|). These representations are nonequivalent, unless - AJ. 6 Z for 
all indices k. 

Given the representation ( |HD of the quantum algebra A in C°°(T m ), any polynomial 
Hamiltonian 7i(Ik) of a completely integrable system is uniquely quantized as a Hermitian 
element Ti(Ik) = 7~t{Ik) of the enveloping algebra A of A. It has the countable spectrum 
(^j. Note that, since Ik are diagonal operators, one can also quantize the Hamiltonians 
which are analytic functions on M. m . 
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3 The classical controllable completely integrable system 

A generic momentum phase space of a Hamiltonian system with time-dependent parameters 
is a composite fiber bundle 

n -> s -> R, 

where II — > £ is a symplectic bundle, E -> R is a parameter bundle whose sections are 
parameter functions, and R is the time axis 12], |19| . Here, we assume that all bundles are 
trivial and, moreover, their trivializations hold fixed. Then the momentum phase space of a 
completely integrable Hamiltonian system on the symplectic annulus W ([!]) is the product 

n = »x5xiy, (18) 

equipped with the coordinates (t,a x ,Ik,4> k ). It is convenient to suppose for a time that 
parameters are dynamic variables. The momentum phase space of such a system is 

II' = lx T*S x W, (19) 

coordinated by (t,<r x ,p\,Ik,(fi k )- 

The dynamics of a time-dependent mechanical system on the momentum phase space 
n' (|19"D is characterized by a Hamiltonian form 

# s = p x da x + I k d<p k - Hx(t, a», Pli , Ij, (jP)dt (20) 



12] , [18| . For any Hamiltonian form H% (|20|), there exists a unique vector field 7# on n' 



such that 

j H \dt = l, jH\dHY; = 0. 

It defines the first order differential Hamilton equation on n'. 
Note that the Hamiltonian Ti^ (p0|) takes the form 



= Px T x + I k (A k + T x Al)+H, (21) 

where 

A o T = dt ® (d t + T x d x + (A fc + T x A k )d k ) 
is a connection on 1 x 5 x T m — > R which is the composition of a connection 

r = rft® (<9 t + r A (t, o<9 A ) (22) 

on the parameter bundle RxS^R and of a connection 

A = dt® (d t + A k (t, a", <fP)d h ) + da x ® (9 A + Aj(t, a", <j?)d k ) (23) 
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on M x S xT m ^ R x S 0, [T§. 

Bearing in mind that a x are parameters, one should choose the Hamiltonian TL^, ( FT] ) 
affine in momenta p\. Furthermore, in order to describe a Hamiltonian system with a fixed 
parameter function a x = £ A (£), one defines the connection T (|22|) such that 

v r e = o, r A (t,^(t)) = ^e A . 

Then the pull-back 

H = £*F E = J fc rf0 fe - (J fc [A fc (t,e,0 J ) + At(t,^,<me A ] +H(t,e,Ij,P))dt 
is a Hamiltonian form onEx H 7 . Let us put 

H = H- I k A k . 
Then the Hamiltonian form 

H = I k d<f) k - H'dt = I k d<f) k - [I k A k x (t, e, ^)d t e + W^dt (24) 

describes a time-dependent perturbed completely integrable Hamiltonian system on M x W . 
The corresponding Hamilton equation reads 

dth = -d k h.{i 3 d t i\ d t <p k = d k H + A k d£ X - (25) 

In order to make the term 

A = I k A k d t e (26) 



in the perturbed Hamiltonian Ti! (|24]) to a control operator, let us assume that the coeffi- 
cients of the connection ( p3| ) are independent of time. Then, in view of the trivialization 



T^), the second term 

da x ^(d x + A k (^,^)d k ) 

of the connection (^) can be seen as a connection on the fiber bundle S x T m —>■ S. Let the 
dynamic Hamiltonian TC be independent of action variables with some I indices a, b, c, . . ., 
and let the perturbation term A ( |26| ) be independent of the action and angle variables with 
the rest indices k, . . .. Then the Hamilton equation (|25|) falls into the two independent 
equations 

d t I k = 0, d t <p k = d k H, (27) 
(a) d t I a = -I b d a A b M\ (b) d t( p a = A a x d t e- (28) 

The first equation (p7|) keeps the form of the Hamilton equation ((^) of an autonomous 
completely integrable system, while the second one is the control equation as follows. 
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Let us rewrite the equation (|28|b) as the countable system of equations 



(29) 



for functions Vv<0 ©• ^ ne hand-side of these equations is a multidimensional Fourier 
series with time-dependent coefficients. Therefore, the equations (|2~9"D for all collections of 
I integers (n a ) make up a countable system of ordinary linear differential equations with 
time-dependent coefficients: 

d^ {na) =iM { ^ a) {e)d t e^n b ). 

Its solution with the initial dates <j) a (0) can be written formally as the time-ordered expo- 
nential 

( n b), 



U(t) = Texp 



M x {r{t'))d t edt' 



Texp 



Mxia^da* 



C([o.*]) 



(30) 



T^, A glance at the evolution operator U(t) fl3~0"D shows that solutions ip(n a )(t) of the 
equations (|29"D are functions of a point a of the curve £ : R — > S 1 in the parameter space S. 
Substituting this solution into the equation (EBk), we obtain the system of I ordinary 



d a A\. 



linear differential equations with time-dependent coefficients: 

dj a = -L b Xa ^ (na) (t),e(md t e, l\ 

Its solution is given by the time-ordered exponential 
I a (t) = U(t) b J b (0), 



U(t) = Texp 



- I L x (^ na) (t'),e(t'))d t edt' 



Texp 



L x (iP (na) (a),cT»)da x 

eao,*]) 

S in the parameter space 



This solution is also a functions of a point a of the curve £ : K. 
S. 

Consequently, the equation (|28|) can be regarded a control equation. In particular, one 
can choose a parameter function such that the trajectory of the perturbed system 
does not leave a compact neighbourhood of the invariant torus T m of the initial completely 
integrable Hamiltonian system. 
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4 The quantum controllable completely integrable system 

Let us quantize the classical controllable completely integrable system on the momentum 
phase space K x W in previous Section. 

The manifold W = R x W, coordinated by (t,I k ,<fi k ), is provided with the Poisson 
structure 

{/, /'} = d k fd k f - d k fd k f, /, /' e C°°(W). 

This is the direct product of the symplectic structure Q (0) on the symplectic annulus W 
(|TJ) and of the zero Poisson structure on the time axis R. In particular, the Poisson algebra 
(C°°(W), {, }) of smooth real functions on W is the Lie algebra over the ring C°°(R) of 
smooth real functions on R. In order to quantize the Poisson manifold (W, {, }), we there- 
fore can essentially simplify the general procedure of instantwise geometric quantization of 
time-dependent Hamiltonian systems in Ref . . Namely, we repeat geometric quantization 
of the symplectic manifold (W, fl) in Section 2, but replace functions on T m with those on 
R x T m . 

Let us choose the angle polarization spanned by the vectors d l . The corresponding 
quantum algebra A C C°°(W') consists of affine functions 

f = a i (t,cfS)I i + b(t,<y) 

of action coordinates /j. It is represented by the first order differential operators 

/ = -i a % - ^ditf - dXi + b, KE R, (31) 

in the space C°°(RxT m ) of smooth complex functions on RxT m . Given different collections 
of real numbers (Aj) and (A-), the representations fl3"Ip are nonequivalent, unless Aj — A^ eZ 
for all indices i. The carrier space C°°(R x T m ) is provided with the structure of the 
pre-Hilbert C°°(R)-module with respect to the nondegenerate C°°(R)-bilinear form 

/ 1 \ m r 

= [^) J W , ^ e C°°(R x T m ). 

The basis for this pre-Hilbert module is made up by the pull-backs onto R x T m of functions 
ip(n r ) (0) on T m . They are the eigenvectors of the action operators 

I k = -id k - A fc , = ( n k - h)^(n r )- 

As in previous Section, let us assume that the dynamic Hamiltonian TC is independent 
of action variables with some / indices a,b,c, . . ., and let the perturbation term A ( f2"6[) be 
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indepependent of the time and the action-angle variables with the rest indices k, . . ., 
i.e., the perturbed Hamiltonian TC' ( JMp takes the form 

W = Kl{e^ b )d t ei a + H{I J ). 

The perturbation term 

of this Hamiltonian is an element of the quantum algebra A, and is quantized by the 
operator 

A = -iA a p d t fd a - l -d a {Kl)d t e - K^e t e. 

The (polynomial or analytic) dynamic Hamiltonian TC(Ij) is quantized as in Section 2, i.e., 
H = H(T 3 ). 

Since the operators A and 7i mutually commute, the corresponding quantum evolution 
operator reduces to the product 



Texp 



-i / n'df 



U 1 (t)oU 2 (t)=Texp 



t 




t 


-i J ndt' 


o Texp 


-i J Adt' 











(32) 



The first operator in this product is the dynamic evolution operator of the quantum com- 
pletely integrable Hamiltonian system. It reads 



Ui(t)ip {nj) = exp[-iH(nj - A,,-)*]^)- 



(33) 



Its eigenvalues are countably degenerate. Recall that the operator (^) acts in the C°°(R)- 
module C°°(RxT m ). Its eigenvalues are smooth complex functions on R, and its eigenspaces 
are C°°(R)-submodules of C°°(R x T m ) of countable rank. 
The second multiplier in the product ( |3!f) is 



U 2 (t) = Texp 
= Texp 



{-A a j<f> b ,^(t'))d a - ^A^exo) +i\ a K^\^{t'))}d t edt' 



{-A;(0 b , <r»)d a - -d a A%(<p b , a») + t\ a A a ^ b , <j»)W 



(34) 



It acts in the eigenspaces of the operator U\(t). For instance, such a space is exemplified 
by the pre-Hilbert C°°(R)-submodule Sq C C°°(R x T m ) whose orthonormal basis is made 
up by functions i/j( na ) for all collections of integers (n a ). Written with respect to this basis, 
the operator [^(t) acts on £ as a matrix of countable rank. 
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A glance on the expression (R3) shows that, in fact, the operator U 2 (t) depends on 



the curve £([0,1]) C S in the parameter space S. One can treat it as an operator of 
parallel displacement with respect to a connection in the C°°(S)-module of smooth complex 
functions on £ x T m along the curve £ || [13|, For instance, if £([0, 1]) is a loop in S, 
the operator Ui (|34]) is the geometric Berry factor. In this case, one can think of U 2 as 
being the holonomic control operator. 

It should be emphasized that the adiabatic assumption has never been involved. 
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